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SUMMARY 

Using high-frequency approximations, the secondary pattern of a reflector 
antenna can be calculated by numerically evaluating a radiation Integral 
I(u,v). In recent years, tremendous effort has been expended to reducing 
I(u,v) to Fourier Integrals. These reduction schemes are Invariably reflector 
geometry dependent. Hence, different analyses/computer software development 
must be carried out for different reflector shapes/boundaries. The purpose of 
this note Is to point out, that, as the computer power Improves, these reduc- 
tion schemes are no longer necessary. Comparable accuracy and computation 
time can be achieved by evaluating I(u,v) by a brute-force FFT described In 
this note. Furthermore, there Is virtually no restriction on the reflector 
geometry by using the brute-force FFT. 


RADIATION INTEGRAL 

In calculating the radiation pattern of a reflector antenna, the central 
step Is to evaluate the following radiation Integral: 

I(u.y) -// J(x.y) ux+vy+wz) dx dy < 

l 

where, 

(e, 4 >) observation direction 
k = 2ir/\ = wave number 

u = sin e cos 4 > 

v = sin e sin «j> 

w = cos e = -yi - (u2 + v 2 ) 




z = f(x,y) (equation describing the reflector surface) 

E = projection of the reflector on the xy plane (fig. 1(b)). 

We note that equation (1) Is not a two-dimensional Fourier Integral because of 
the presence of the factor exp(jkwz), which reads explicitly 


exp(jkVl - (u 2 + v 2 ) f(x,y)) (2) 

In the past years, there has been a tremendous amount of effort to convert 
the radiation Integral (1) to a Fourier Integral. A popular method Is the 
so-called "p-serles expansion," In which the Integral In equation (1) Is 
expressed as 


I(u,v) = ^ (cos e - cos e Q ) p I p (u,v) (3) 

Here eg Is the polar angle of the main beam (fig. 1(a)). Each Integral 
Ip(u,v) In equation (3) is Indeed a Fourier Integral. Several clever schemes 
can be used for Its evaluation (refs. 1 to 6). As an alternative to the 
p-serles expansion, we may first calculate the aperture field over a planar 
surface by GO or GTD (refs. 7 and 8). Over the planar surface, z Is equal to 
a constant and the Integral In equation (1) Is reduced to a Fourier Integral. 
Needless to say, either the p-serles expansion or GO/GTD involves sophisticated 
mathematical manipulations. 

While the above efforts are gratifying, they all suffer a drawback, 
namely, the analyses are heavily geometry dependent. Usually a given analysis 
falls If the reflector surfaces Is changed from, for example, parabolic to 
hyperbolic, or the reflector boundary from circular to trapezoidal. As a con- 
sequence, one needs to develop different versions of computer code for differ- 
ent reflector geometries. 


BRUTE-FORCE FFT 


The purpose of this note Is to point out that, because of the rapid pro- 
gress of computer power. It Is no longer advantageous to reduce equation (1) 
to a Fourier Integral before going to a computer. We evaluate equation (1) by 
a brute-force FFT In the following manner: 


(a) Represent I(u,v) In terms of a sine series, namely, 


I(u,v) 



t s1n(ir( uA - m) 
A mn ir( uA - m) 


s1n(ir(vB - n)) 
ir(vB - n) 


(4) 


2 



where 


t M_ 
mn = 4MN 


j2ir(ux. ) +vy. ] ) 


p e J2ir(uA-m+vB-n) 
mn 


F mn = ^ our ^ er ser ^ es coefficients of J(x,y) 


e kjwz 


(5) 


Note that the coefficients {I mn } are values of I(u,v) at sampling points 

u = u m = m/A, v = v n = n/B (6) 

(b) Use the FFT to find {I mn } from equation (1). Since equation (1) 

Is not a Fourier Integral, we must apply the FFT once for each and every I mn . 
Thus, we apply the FFT 2M by 2N times Instead of once, If equation (1) was a 
Fourier Integral. Fortunately, modern computers can perform the FFT operation 
very efficiently. 


ACCURACY AND COMPUTATION TIME 

To Illustrate the accuracy and efficiency of the brute-force FFT method, 
the secondary pattern of an offset parabolic reflector with a circular aperture 
(fig. 2) was computed using the brute-force FFT, GO (ref. 7), GTD (ref. 7), 
and Fourler-Bessel (ref. 1) methods. In each method, a 48 by 48 FFT grid 
(fig. 1(b)) was taken and the secondary pattern was calculated at 101 observa- 
tion points. As expected, the secondary patterns calculated by each method 
are In good agreement. Taking a closer look, the gain and peak sldelobe levels 
for the first six sldelobes are tabulated for each method In table I. These 
values are In good agreement with the results calculated using the Jacobl- 
Bessel series method (ref. 5). Another example with an off-focus feed Is shown 
In figure 3, where good agreements among all methods are again observed. 

Table II shows the computation times of the methods on a CYBER 175 com- 
puter. We note that the computation time of the brute-force FFT Is of the 
same order of magnitude as the computation times using GO, GTD, and the 
Fourler-Bessel series method. For nonparabollc reflectors, the coefficients 
of the p-serles method must be computed for every observation point. Hence, 
the computation times using the p-serles are of the same order as the brute- 
force FFT. The computation times using GO and GTD may also Increase due to 
the possibility of more specular points present for a nonparabollc reflector. 


CONCLUSION 

1. Unlike other methods described In the literature, the brute-force FFT 
Is most general In the sense that It applies to reflector with arbitrary shape 
and boundary. It Is an Ideal method for developing a user-friendly general 
purpose reflector computer program. 

2. Within the high-frequency approximation used In all reflector 
analyses, the brute-force FFT Is just as accurate as the other methods. 
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3. Its computation time Is In the same order of magnitude as those of 
other methods. In the worst case (a parabolic reflector with circular aper- 
ture), It Is three times slower than other methods that take advantage of the 
special geometry and use clever mathematical manipulations. 

4. The above time comparison is based on the fact that the FFT Is per- 
formed using a software subroutine. For the Cyber 175 computer, the FFT com- 
puted using an array processor Is five times faster than that using a software 
subroutine. Our message Is that, as computer power Improves, the brute-force 
FFT for reflector analysis will become more and more attractive. 
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TABLE I. - GAIN AND SIDELOBE LEVELS OF OFFSET PARABOLIC 


REFLECTOR COMPUTED USING SEVERAL DIFFERENT METHODS 



Jacobl-Bessel , 

FFT 

GTD 

GO 

Fourler-Bessel , 


ref. 5 


ref. 7 

ref. 7 

ref. 7 

Gain (dB) 

48.28 

48.25 

48.31 

48.29 

48.22 

1st SL 

28.70 

28.94 

28.48 

28.49 

27.98 

2nd SL 

22.50 

22.48 

22.64 

22.78 

21.92 

3rd SL 

18.32 

18.85 

17.72 

17.63 

17.76 

4th SL 

15.21 

15.55 

15.80 

14.93 

14.64 

5th SL 

12.67 

12.76 

12.09 

8.54 

12.14 

6th SL 

10.55 

11.15 

11.30 

9.96 

10.01 


TABLE II. - EXECUTIVE TIME OF SEVERAL 
DIFFERENT METHODS ON A CYBER 175 COMPUTER TO 
COMPUTE SECONDARY PATTERN 


























(a) Main beam reference (polar angle). 



(b) Fourier aperture grid. 
Figure 1. 
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Figure 2. - Offset parabolic reflector configuration. 



Figure 3. - Far-field antenna pattern comparison. 
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